MULTIPLICITIES OF THE MOST SINGULAR POINT ON 
SCHUBERT VARIETIES ON GL(N)/B FOR n = 5,6 
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Abstract. We calculate using Macaulay 2 the multiplicities of the most sin- 
gular point on Schubert varieties on Gl(n)/B for n = 5, 6. The method of 
computation is described and tables of the results are included. 



1. Introduction 

In this paper, we compute, using Macaulay 2, the multiplicity of the most singu- 
lar point on Schubert varieties of the flag manifold Gl(n)/B for n — 5 and n = 6, 
using a description of the preimage of Schubert varieties in GL(n) first given by 
Fulton j3| and more recently developed by Knutson and Miller [Sj. Results of 
Krattenthaler 6 , Rosenthal and Zelevinsky j^j, and Lakshmibai and Weyman [H] 
give combinatorial and determinantal formulas for multiplicities (at all points) of 
Schubert varieties on the Grassmanian. Furthermore, the singular loci of Schubert 
varieties on flag manifolds have been much studied, with known results collected 
in PP, starting with the fundamental result of Lakshmibai and Sandhya 7 that a 
Schubert variety indexed by the permutation w is singular iff w contains either the 
pattern 1324 or the pattern 2143. However, no results for the multiplicities of the 
Schubert varieties appear to be known on the full flag variety. 

In the next section, we briefly define the objects under study and outline some 
basic results about them. A detailed introduction to Schubert varieties can be 
found, for example, in part III]. This section also serves to fix our conventions; 
the differing choices made by different authors in the subject can cause significant 
confusion. In particular, our convention for indexing Schubert varieties is opposite 
to the convention used in p^, which now seems to be fairly standard. The third 
section describes our algorithm, and the fourth section demonstrates this algorithm 
for w — 1324. Two appendices give our code and the results of our computations. 

2. Definitions and Conventions 

A (complete) flag T in C" is a sequence of subspaces {0} = F C F\ C F 2 C 
• • • C F n = C" such that the subspace Fi has dimension i. Fixing a basis for C n , 
we can represent T by a matrix M as follows. For each component Fi of J 7 , pick a 
vector rrij e Fj \ Fj_ l5 and write m, as the i-th row of M. Note that Fi will be the 
span of the first i rows of M; in particular, M is in invertible since F n = C". 

This representation of the flag is clearly not unique, as it involves repeated choices 
of vectors. To be precise, both multiplying any row of M by a constant and adding 
any row to a subsequent row of M leaves the flag unchanged. This is equivalent 
to multiplying M on the left by a lower triangular matrix. We can now give the 
structure of an algebraic variety to the set of flags; namely, it is the quotient of 
G = Gl(n) by the group of lower triangular matrices B~ acting on the left, and we 
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denote this variety by B \G and call it the flag variety. We also have a natural 
map 



sending a matrix to the flag it represents. 

We can, however, pick a standard representation N(T) for each flag T as follows. 
Pick a matrix M that represents T. Now take the leftmost nonzero entry in the 
first row of M = {rriij}, which we call mu, and use elements of B~ to change M so 
that mu = for k > 1. We then take the leftmost nonzero entry of the second row, 
place Os in all entries below it, and repeat for all rows in order to finally get N(T). 
Unfortunately, we cannot identify B~\G with these representations of flags except 
as a set, since this process of row reduction destroys the geometry of B~\G. 

Let the group B + of upper triangular matrices act on G by right multiplication; 
by associativity of multiplication, this commutes with the left action of B~ and 
therefore gives an action on B~\G. This action can be thought of as adding any 
column of a matrix to a column to its right. Under these two actions, any matrix 
can be sent to a unique permutation matrix W, so we can index orbits of the right 
action of B + on B~\G by permutations. In particular, a flag T is in the orbit of 
tt(W), where W is the permutation matrix with Is in the leftmost nonzero entries 
of each row in N(T). Let X w denote the orbit of the flag 7r(VU -1 ), where W is the 
permutation w written as a matrix. X w is known as a Schubert cell. The process 
of choosing the representative N(!F) for each flag does preserve geometry locally 

on X w , so X w is isomorphic to a( 2 ^~ 1(w \ where l(w) is the length of a shortest 
expression of w as a product of adjacent transpositions, or, equivalently, the number 
of inversions in w. 

We denote by Y w the closure of X w in B~\G; it is known as a Schubert variety. 
For permutations v, w S S n , let v y w if l(v) > l(w) and v = tw for some transpo- 
sition t. The transitive closure of the relation >- is known as the Bruhat order, for 
the remainder of this paper, v > w for v, w G S n means that v is greater than w in 
this partial order. It is a classical result that 



Note that the unique O-dimensional cell X WQ , where u>o — n- ■ -21 is contained 
in every Schubert variety. 

Given a variety X, the multiplicity at a point p of J, which we will denote 
multp X is the degree of the projective tangent cone Proj(gr mp Ox.p) as a subvariety 
of the projective tangent space Proj(Sym* m p /m^), or, equivalently, if the Hilbert- 

Samuels polynomial of Ox, P is written a n x n +a„_ix™~ 1 H ha , mult p X = n\a n . 

For a nonsingular point p, mult p X = 1; at a singular point q, mult 9 X > 1 and 
measures roughly how singular X is at q. Slightly more precisely, the multiplicity 
counts how many times a generic hyperplane cuts through X in a neighborhood of 

Since a Schubert variety is invariant under the right action of B + , its multiplicity 
must remain constant on B + orbits, or Schubert cells. Moreover, since for v >- 
u>, there exists a P 1 with one point in X v and the remaining points in X w , by 
scmicontinuity, multiplicities must be nondecreasing with respect to Bruhat order. 
X Wo must therefore be the most singular point of Y w and the multiplicity of Y w 
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there measures how singular Y w gets. In particular, a multiplicity of 1 at X-, 
indicates that Y w is smooth. 



Since multiplicity is a local property, we can calculate it after restricting to an 
affine neighborhood of X Wn in B~\G. A natural candidate is Sl Wo , the orbit of 
Wo (defined to be wq written as a permutation matrix and considered as a flag) 
under the right action of B~ . (In general, Q, w is defined as the orbit of the flag 
iriW -1 ) under the right action of B~~ (rather than B + ) and is known as a dual 
Schubert cell.) Locally on £l Wo , the map n : G — > B~\G has a section a, namely 
the map that sends a flag T to N(J-). This identifies Cl wo with the matrices with 
Is on the main antidiagonal and Os to the right and below; X wo is mapped to the 
permutation matrix Wq- Since a is a local section, Y w is locally isomorphic in a 
neighborhood of X Wo to tt" 1 (Y w ) n a(Q Wo ). 

Now we need to find equations defining n (Y w ). Fix a permutation w G S n . 
Let R(w) = [rij(w)] be the integer matrix with Vij(w) — #{w (k) < i,k < j}. 
For any matrix M, let denote the submatrix consisting of the first i rows and 
first j columns. Then, for any invertible matrix M with 7r(M) G X w , the rank of 
the submatrix My will be r»j. The proof of this claim is as follows. Note that the 
rank of Mij for any i and j does not change under multiplication by B~ on the 
left, since the effect of multiplication by b G B~ on the first i rows is the same as 
that of multiplying by an element of Gl(i), namely the submatrix of b consisting 
of the first i rows and columns. Therefore, the claim can be verified on matrices of 
the form N^J 7 ) for T G X w , where it is trivial. 

It is a nontrivial combinatorial fact that Bruhat order can be equivalently defined 
by v > w if Tij(v) < rij{w) for all Therefore ir (Y w ) consists of all invertible 
matrices M satisfying the rank conditions rk(My) < rij(w). Now let Z = [zij] be 
a matrix of indeterminates, and I w be the ideal generated by all size 1 + nj(w) 
minors of Z^, for all i and j. (As shown in j^j, where I w was originally defined, a 
small subset of the minors suffices to generate I w , but we will not use this fact.) By 
the above statement, it is clear that I w vanishes precisely on the points of tt~ 1 (Y w ); 
Knutson and Miller [H] prove that I w is in fact radical, so I w = /(7r _1 (Y w )). Now, 
by sending siij to for entries below the main antidiagonal and 1 for entries on the 
main antidiagonal, we have the ideal J w for tt^ 1 (Y w ) n a(Cl Wo ) as a subvariety of 



In our coordinates for a(Q Wo ), X W[) corresponds to the point where = for all 
i and j, or equivalently the maximal graded ideal m = (zij). The projective tangent 
cone of Y w is therefore Proj(gr m S/J w ), where S = For a polynomial / G S, 

let fd be the homogeneous part of degree d, and let s(f) be the smallest number 
such that / s(/) ^ 0. Let J' w = (/ s (/)|/ G J w ); then gr m S/J w = S/J' w , since any 
polynomial / is sent to / (mod m s+1 ), of which f s is a representative. Degree is 
invariant under Grobner deformation and easiest to calculate on monomial ideals, 
so, to calculate the multiplicty of Y w , or, equivalently, the degree of S/J' w , it suffices 
to calculate the degree of S/ in( J' w ) under any (graded) term order. This is the same 
as taking the initial ideal of J w , taking care to use a term order which places the 
lowest degree term first. 

In practice, on Macaulay 2, one does this by homogenizing the generators of J w 
using a new variable t (or by sending z^j to t rather than 1 for entries on the main 



3. Explanation of the Algorithm 
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diagonal in passing from I w to J w ), and using a term order that refines the partial 
order by degree in t. (See Prop 15.28] for a proof that this is equivalent.) We 
can then compute the initial ideal, send t to 1, and calculate the degree. 



4. Example 

As an example, we calculate the multiplicity of Y w at X WQ for w — 2143, the 
smallest nontrivial example. Here, we have 



W 



1 

1 





Therefore, we have the rank matrix 








Exactly two of the rank entries give minors in I u 
z\\ G I w , and r 33 = 2, which gives 



namely m = 0, which gives 



Zll Z U Z W 

det z 2 i z 2 2 z 2 3 = 0. 
_z 3 i z 32 z 33 

Setting z 33 = and z 32 = z 23 = 1, we have J w = (zn, —z\\ + z 12 z 31 + z 21 z X3 - 
z 3 \z 22 z 13 ) = (z n , z\ 2 z 31 + z 21 zi 3 - z 31 z 22 z 13 ) . Then J' w = {zn,z 12 z 31 +z 21 zi 3 ), and 
it is clear that the multiplicity of Y w at X Wo is 2. For purposes of illustration, we 
carry out the remainder of the algorithm. Homogenizing the generators of J w gives 
us (z\\,tz\ 2 z 3 \ + tz 2 \Z\ 3 — z 3 \z 22 zi 3 ), and one possible appropriate initial ideal is 
{z\i, tz 2 \Z\ 3 ). Sending t to 1, we get (zn, z 2 \Z\ 3 ) , which has degree 2; therefore, we 
conclude that multj^ , Y w — 2. 
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Appendix A. Maple and Macaulay 2 code 

The following is our Maple code for generating the rank matrices. 

with(combinat) ; 
with(linalg) ; 

interface (prettyprint=f alse) ; 
n:=5; 



sum (matrix (n , n , 

(i, j)->'if ' ((j>=k) and (i>=op(k, pPerm)), 1, 0)), k=l..n); 
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The Maple output was sent to a file, and all square brackets [] were converted 
to curly brackets {} for Macaulay 2. Macaulay 2 code is as follows. Note that the 
first five lines must be changed for each n. 

R=QQ[t,xll,xl2,x21,xl3,x22,x31,xl4,x23,x32,x41, 
MonomialOrder=>Eliminate 1] ; 

G = matrix({{xll,xl2,xl3,xl4,t}, 
{x21,x22,x23,t,0>, 
{x31,x32,t,0,0}, 
{x41,t, 0,0,0}, 
{t, 0,0, 0,0}}) 

S=QQ [Xll,xl2,x21,xl3,x22,x31,xl4,x23,x32,x41] ; 

f =map(S,R,{l ,xll,xl2,x21,xl3,x22,x31,xl4,x23,x32,x41}) ; 

n=5; 

# Mlist = (paste from Maple output) 

# compute J_w 
Hist = apply (Mlist, 

M->trim(sum(f latten(f or i from to n-1 list 

for j from to n-1 list 
minors(M_(i, j)+l, 

submatrix(G, {0..i}, {0. . j}))))) ; 

GBlist = applydlist, gb) ; 

# gives in(J_w) 

LTlist = apply(GBlist, GB -> leadTerm(gens (GB) ) ) ; 

# gives in( J~\prime_w) 

ELTlist = apply (LTlist, LT->f (LT) ) ; 

# outputs degrees 

Dlist = apply (ELTlist, LT -> degree (ideal (LT) ) 

The output of the last line gives the degrees. Output from other lines, such as 
the initial ideals, could also be of interest. These computations run quite quickly; 
in fact, generating the ideals J w was by far the slowest step, and the Grobner basis 
computation took only a few seconds for n = 6. 
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Appendix B. Computational Results for n — 5 and n = 6 

We have listed the permutations by multiplicity of Y w at X WQ . For each multi- 
plicity, permutations are listed in lexicographic order. 



Multiplicity 


Permutations 


5 


14325 


















3 


13425, 


14235, 


24153, 


31524 












2 


12435, 


13245, 


13254, 


13524, 


14253, 14352, 


15324, 


21354, 


21435, 




21453, 


21534, 


21543, 


23154, 


24135, 


24315, 


24351, 


25143, 


31254, 




31425, 


31542, 


32154, 


32514, 


32541, 


41325, 


42153, 


51324, 


52143 


1 


12345, 


12354, 


12453, 


12534, 


12543, 


13452, 


13542, 


14523, 


14532, 




15234, 


15243, 


15342, 


15423, 


15432, 


21345, 


23145, 


23415, 


23451, 




23514, 


23541, 


24513, 


24531, 


25134, 


25314, 


25341, 


25413, 


25431, 




31245, 


31452, 


32145, 


32415, 


32451, 


34125, 


34152, 


34215, 


34251, 




34512, 


34521, 


35124, 


35142, 


35214, 


35241, 


35412, 


35421, 


41235, 




41253, 


41352, 


41523, 


41532, 


42135, 


42315, 


42351, 


42513, 


42531, 




43125, 


43152, 


43215, 


43251, 


43512, 


43521, 


45123, 


45132, 


45213, 




45231, 


45312, 


45321, 


51234, 


51243, 


51342, 


51423, 


51432, 


52134, 




52314, 


52341, 


52413, 


52431, 


53124, 


53142, 


53214, 


53241, 


53412, 




53421, 


54123, 


54132, 


54213, 


54231, 


54312, 


54321 







For n = 6, we have omitted the permutations of multiplicities 1 and 2 for reasons 
of space. (There are 366 nonsingular Schubert varieties and 207 with multiplicity 
2 at X Wo .) These can be inferred from the remainder with the help of the pattern 
avoidance criterion. 



Multiplicity 


Permutations 


14 


154326 














10 


153426 














9 

8 


145326, 
321654 


154236 












7 


135426, 


143526, 


152436, 


153246, 


254163, 416325 






6 


145236, 


132546, 


214365 










5 


125436, 


135246, 


142536, 


143256, 


143265, 143625, 


146325, 


153264, 




154263, 


154362, 


164325, 


215436, 


251364, 251436, 


251463, 


253164, 




254136, 


254316, 


254361, 


314625, 


315426, 316425, 


413625, 


415326, 




426153, 


514326, 


614325 










4 


153624, 


152346, 


134526, 


214635, 


215364, 215463, 


216435, 


231564, 




231654, 


241365, 


243165, 


245163, 


312645, 312654, 


314265, 


321564, 




321645, 


326154, 


351426, 


351624, 


413265, 416235, 


421653 




3 


124356, 


124365, 


124536, 


125346, 


132564, 132645, 


132654, 


134256, 




134265, 


134625, 


135264, 


136425, 


142356, 142365, 


142635, 


145263, 




145362, 


146235, 


152364, 


152463, 


153462, 163425, 


164235, 


214356, 




214536, 


215346, 


216453, 


216543, 


231546, 235164, 


241536, 


241563, 




241635, 


241653, 


245136, 


245316, 


245361, 246153, 


251346, 


251634, 




251643, 


253146, 


253416, 


253461, 


264153, 312546, 


314526, 


315246, 




315264, 


315624, 


316245, 


316254, 


316524, 321546, 


325164, 


341625, 




351264, 


351642, 


352164, 


352614, 


352641, 361524, 


425163, 


412635, 




413526, 


415236, 


416253, 


416352, 


421635, 423165, 


426135, 


426315, 




426351, 


431625, 


432165, 


513426, 


514236, 524163, 


531624, 


613425, 




614235, 


624153, 


631524 
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